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Recent studies have shown that gamma-ray pulsar Ught curves are very sensitive to the geometry of the pulsar 
magnetic field. Pulsar magnetic field geometries, such as the retarded vacuum dipole and force-free magneto- 
spheres have distorted polar caps that are offset from the magnetic axis in the direction opposite to rotation. 
Since this effect is due to the sweepback of field lines near the light cylinder, offset polar caps are a generic 
property of pulsar magnetospheres and their effects should be included in gamma-ray pulsar light curve mod- 
eling. In slot gap models (having two-pole caustic geometry), the offset polar caps cause a strong azimuthal 
asymmetry of the particle acceleration around the magnetic axis. We have studied the effect of the offset polar 
caps in both retarded vacuum dipole and force-free geometry on the model high-energy pulse profiles. We find 
that, compared to the profiles derived from symmetric caps, the flux in the pulse peaks, which are caustics 
formed along the trailing magnetic field lines, increases significantly relative to the off-peak emission, formed 
along leading field lines. The enhanced contrast produces improved slot gap model fits to Fermi pulsar light 
curves like Vela, with vacuum dipole fits being more favorable. 
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The Fermi Gamma-Ray Space Telescope has had 
a major impact on pulsar physics with the discov- 
' ery of over 100 gamma-ray pulsars comprising three 
^ populations: young radio-loud pulsars, young radio- 
quiet pulsars and millisecond pulsars f22l. The wide 
variety of light curve types, viewing geometry and 
spin-down luminosity encompassed by this large group 
of sources provides an unprecedented opportunity to 
explore and constrain pulsar emission and magnetic 
field geometry. There have also been significant ad- 
vances recently in numerical simulation of pulsar mag- 
netosphere models [53] U [TS] [T5] that define the high- 
altitude magnetic field structure critical to the emis- 
sion gap models. Since the radiation in the outer mag- 
netosphere high-energy emission models, such as outer 
gap ^] and slot gap TS], occurs along the bound- 
ary of the open field region, the predicted gamma-ray 
light curves are sensitive to the magnetic field struc- 
ture [H]. In particular, the sweepback of the magnetic 
field lines near the light cylinder, due to retardation 
and currents, causes an offset of the polar cap in the 
direction opposite to the rotation. In the slot gap 
(SG) model, the polar cap offset produces two main 
effects on the emission and the resulting light curves: 
an azimuthal asymmetry in the SG emission, and a 
change in the phase lag between the first gamma-ray 
peak and the radio peak. We present a study of SG 
model light curves for the two extreme cases of vac- 
uum dipole and force-free pulsar magnetospheres to 
examine such effects. 



2. PULSAR MAGNETOSPHERE MODELS 

The global structure of a pulsar magnetosphere is 
a presently unsolved problem. An analytic expression 
exists only for the structure of the magnetosphere of 
the vacuum rotator (the Deutsch solution, [7]), but 
real pulsars are certainly not in vacuum since electrons 
and positrons are copiously produced due to the high 
surface electric fields induced by rotation. The other 
extreme, known as force-free models [S][23][51] that 
specify the configuration of fields, charge and currents 
that solve Maxwell's Equations in the approximation 
of ideal-MHD (E • B = 0) neglecting plasma inertia, 
is also well determined. Models with finite conductiv- 
ity that bridge these two extremes have been studied 
very recently [TS][Tn]. These simulations show that 
there is a range of magnetic field structures, current 
distributions and spin-down power that lie between 
the vacuum dipole, with maximum accelerating elec- 
tric field but no charges, and the force-free solutions, 
with charges and current but no acceleration. 

One property that all these solutions have in com- 
mon is sweepback of the magnetic field lines near the 
light cylinder opposite to the direction of rotation, 
as the poloidal dipole field transitions into the pulsar 
wind (or electromagnetic wave in the vacuum case). 
The magnetospheres of all these solutions divide into 
open field lines, those which cross the light cylinder, 
and closed field lines, those that close within the light 
cylinder. Tracing the bundle of open field lines to 
their footpoints on the neutron star surface defines 
the polar cap (PC). Due to the sweepback, the PCs 
are asymmetric with their centers shifted backward 
in phase relative to the dipole axis [11110] ■ Examples 
of such offset PCs are shown in Figure 1 for the re- 
tarded vacuum and force-free models. One can see 
that the PC rims in both models are shifted toward 
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Figure 1: Shape of polar caps for retarded vacuum dipole (blue) and force- free (red) magnetospheres for different 
pulsar inclination angles a. 



the trailing side, the amount of shift being a function 
of incUnation angle, with the force-free model PCs 
shifted significantly more than the those of the vac- 
uum model (see also [3|). In addition, the force- free 
PCs are larger than those of the vacuum, especially 
for small inclination, because the volume of open field 
is larger. The notches that appear in vacuum PC rims 
for the higher inclination angles largely disappear in 
the force-free model PCs. 



3. MODEL LIGHT CURVES FOR OFFSET 
POLAR CAPS 



To explore how the magnetic field structure and off- 
set PCs influence gamma-ray pulsar light curves, we 
have generated model light curves using a geometrical 
version of the SG, often referred to as the two-pole 
caustic (TPC) model [5|. The physical SG has its ori- 
gin in polar cap pair cascades that screen the acceler- 
ating parallel electric field above a pair formation 
front (PFF) [2] [II]- Approaching the last open field 
line, assumed to be a perfectly conducting boundary 
where iSy vanishes, the PFF moves to higher altitude 
as particles must accelerate over a longer distance to 
radiate pair-producing photons. A narrow gap of un- 
screened -Eji forms along this boundary, encompassing 
field lines along which the potential drop necessary for 
pair production is never achieved. The electrons in the 
SG keep accelerating to high altitude, radiating cur- 
vature and possibly synchrotron and inverse Compton 
emission [T51 [T^ . 



3.1. Asymmetric Slot Gap Emission 

The original SG model -Ey [2j had azimuthal asym- 
metry because field lines curved toward the rota- 
tion axis rj ('favorably curved') accelerate one sign 
of charge (electrons in the $7 • B > case) and field 
lines curved away from the rotation axis ('unfavorably 
curved') accelerate the other sign of charge. Inclu- 
sion of frame-dragging |17| introduced a larger favor- 
able contribution that significantly reduced this 
azimuthal asymmetry. Both of these -Ey solutions 
assumed centered dipole PCs. Recently Harding & 
Muslimov [l3] [Mj presented a solution for in non- 
dipolar magnetic fields having PCs with offsets of ar- 
bitrary degree and direction. 
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where Bq is the surface magnetic field strength at the 
magnetic pole, rj = r/R is the dimensionless radial 
coordinate in units of stellar radius, R, a ^ e cos((/) — 
00 ) is the parameter characterizing the distortion of 
polar field lines, and is the magnetic azimuthal 
angle defining the meridional plane of the offset PC. 
In this magnetic field, the PC angle has an azimuthal 
dependence. 
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respectively, where x = r/i?LC is the radial distance 
in units of the light cylinder radius, i?LC = c/ft. Since 
the £'|[ has a dependence on the PC angle, its value is 
significantly larger on the offset side of the PC, pro- 
ducing a strong azimuthal asymmetry. 
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Figure 2: Simulated gamma-ray light curves for different combinations of pulsar inclination angle a and viewing angle 
^ with respect to the rotation axis, in a retarded vacuum dipole magnetosphere for symmetric (black) and asymmetric 
(blue) slot gap geometry. The assumed slot gap width is 0.05 of the open field lines and the maximum emission radius 
is 1.2i?LC. The phase of closest approach to the magnetic pole is at 0. 



If we assume that is along the z-axis and that the 
magnetic axis lies in the x-z plane, then the PC offsets 
due to magnetic field sweepback are in the y-z plane 
(00 = '^/'^), so that a = esin0. In this case, the i^y 
of the SG will have an approximate dependence [13] 
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where i?y (0) is the value for the case of a pure dipole 
field. We can quantify the offset parameter e in dif- 
ferent magnetosphere models by the relation 
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where i?tr is the radius of the trailing side of the PC 
and Rpc R{flR/cy^'^ is the radius of a centered 
PC. We find that for the retarded vacuum dipole, e ^ 
0.05 — 0.1 and for the force-free model, e ^ 0.1 — 
0.2, with minimum and maximum corresponding to 
inclination angles of 0° and 90° respectively. 



3.2. Light Curves in Vacuum and Ideal 
Force-Free Magnetospheres 

To simulate SG light curves in magnetospheres with 
offset PCs, we need to know how the asymmetry in 
is related to a corresponding asymmetry in the photon 
emission rate. The Lorentz factors of electrons accel- 
erating in the SG are limited by curvature radiation 
reaction [18] and reach a steady-state 
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where pc is the local field line radius of curvature. If 
the main emission from electrons accelerating in the 
SG conies from curvature radiation, then the photon 
emission rate is equal to the electron energy loss rate 
of 7 oc c>c ii'ii . Thus we can assume that the 

asymmetry in the from Eqn (j3| produces a simi- 
larly asymmetric emission pattern of SG radiation. 

Figure 2 shows an atlas of SG model light curves for 
centered and offset PCs of the vacuum dipole, and dif- 
ferent values of pulsar inclination a and viewing angle 

with respect to the rotation axis. The light curve 
calculations follow that of [3], except that here we 
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Figure 3: Comparison of simulated gamma-ray light curves in retarded vacuum dipole (blue) and force-free (red) 
magnetospheres for asymmetric slot gap geometry. The model parameters are the same as in Figure 2. 



modulate the emission in azimuth according to Equa- 
tion ([3|. For these calculations we assume a slot gap 
width of w = 0.05 and uniform emission along field 
lines in the corotating frame from the neutron star 
surface to maximum radius of r™^'^ = 1.2i?LC: limited 
by a maximum cylindrical radius of r™^'^ = 0.95i?LC- 
The combined phase shifts from aberration and time- 
of- flight, of photons radiated tangent to a magnetic 
dipole field from the polar cap to the light cylin- 
der, nearly cancel those due to field line curvature 
on the trailing edge of the open-field region. Radia- 
tion along such trailing field lines bunches in phase, 
forming sharp emission peaks or caustics in the pulse 
profile. For large a, the caustics from both magnetic 
poles are visible and form two sharp peaks in the light 
curves, whose separation depends on viewing angle 
(. On the leading side, these phase shifts add up 
to spread photons emitted at various altitudes over a 
large range of phases, forming the bridge and off-peak 
emission. From Figure 2, one can see that the light 
curves for centered PCs have high levels of off-peak 
emission, higher than that seen in most Fermi pulsar 
light curves. However, the light curves for asymmet- 
ric SG emission resulting from offset PCs have a much 
lower level of off-peak emission, since emission is mini- 
mum on the leading edge and maximum on the trailing 
edge of the open field volume. As a result, the peak 



level and the peak-to-off-peak ratio is enhanced, but 
does not change the phase lag of the first peak with 
phase (left-hand plot boundary). 

Figure 3 shows the a — C atlas comparing asymmet- 
ric SG light curves for offset PC in the vacuum and 
force-free models for the same model parameters as 
Figure 2. The light curve morphology in the force- free 
magnetosphere is similar to that of the vacuum dipole, 
but the peaks are shifted in phase. The first peak 
of the force-free light curves occurs at later phases 
relative to the magnetic pole crossing, 0.05 - 0.15 de- 
pending on inclination angle, compared to the vacuum 
light curves. Force-free model light curves, computed 
by ^3J for the classic TPC model (rj^t'^ = 0.75i?Lc), 
and by |3] for outer magnetosphere curvature radia- 
tion, show similar delayed phase lags of the first peaks. 
The larger phase lags of the force-free model gamma- 
ray peaks result from two effects: the PCs are larger 
than the vacuum PCs, and the field lines are more 
swept back than in the vacuum dipole magnetosphere. 
Both effects cause a larger offset toward the trailing 
side of the PC, delaying the caustics forming the main 
peaks to later phase (see Figure 1). The radio peak, 
thought to come from cone beam emission centered 
on the magnetic pole at altitudes of several hundred 
km above the neutron star surface, would likely occur 
at earlier phase than that of the magnetic pole due 
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Figure 4: Observed (black) and best-fit model (red) light curves for the Vela pulsar for retarded vacuum dipole (top 
panels) and force-free (bottom panels) magnetospheres. The radio peak is at phase and the vertical blue dotted lines 
mark the phase of the magnetic pole crossing. 



to aberration and retardation, so the total phase lag 
of the first gamma-ray peak to the radio peak will be 
even larger. 

3.3. Test Case: the Vela Pulsar 

The different shapes and radio phase lags of the 
gamma-ray light curves for the vacuum dipole and 
force-free models provide a potentially useful diagnos- 
tic for pulsar magnetosphere structure. As an illus- 
tration, we have fit the Vela pulsar light curve with 
our vacuum dipole and force-free model light curves 
for both symmetric and asymmetric SGs (sec also 6 ). 
We use 30 months of sky survey data, choosing diffuse 
class photons within max[1.6 — 3 log]^Q(i?Gev)7 1-3] de- 
grees of the Vela radio position between 0.1 and 300 
GeV. The light curve is obtained by phase-folding, 
using the Fermi ephemeris, and then displayed in 140 
phase bins of 3000 counts each. We use a Markov 
Chain Monte Carlo maximum likelihood technique 
|25j to determine the best fit a, C, gap width w, max- 
imum emission radius r™^'^ and magnetic pole phase 



lag A^. Our resolution is 1° in a and C for the vacuum 
case and 15° in a (each a separate numerical solution) 
and 1° in ( for the force-free case. The resolution in 
gap width < ui < 0.3 is 0.01 in fraction of the open 
volume, and resolution in maximum emission radius 
0.7i?LC < < 2.0i?Lc is O.IElc- 

The fit results are shown in Figure 4, with the best 
fit model light curves, those with the lowest reduced 
xf, — in red and the observed Fermi light 

curves in back. The best-fit (a, C, w, r"^^^^ A(j), xt) 
values are (52°, 71.5°, 0, 1.0, 0°, 1964) and (64°, 65.5°, 
0.1, 1.2, 18°, 1169) for symmetric and asymmetric vac- 
uum dipole and (30°, 67.5°, 0, 1.5, 332°, 1145) and 
(30°, 64° 0, 1.2, 334°, 1515) for force-free symmetric 
and asymmetric models. In the vacuum dipole model, 
the best fit (a, C) a-re different for the symmetric and 
asymmetric cases since their different levels of off-peak 
emission produce a maximum likelihood in an alter- 
nate location in parameter space. In this case, the 
asymmetric SG model is significantly favored and has 
a smaller P = Q — a and A0, allowing for some aber- 
ration of the radio peak position, which adds to the 
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radio phase lag. The C = 65° in this case also agrees 
with the viewing angle inferred from the X-ray torus 
fit [IH]. In the force-free model, the best fit {a, (, w, 
A(f)) are the same for the symmetric and asymmet- 
ric cases, but the r"^^^ is smaller for the asymmetric 
case. In both force-free fits, the magnetic pole phase 
lags are much larger than those of the vacuum dipole 
models, requiring the radio peak, at phase in the 
plot, to arrive at phase 0.08 later than the magnetic 
pole phase. This would not be expected for a stan- 
dard cone beam unless there was emission from only 
the trailing part of the cone. Additionally, (3 = 37.5° 
in the force-free best-fit models which is much larger 
than the expected width of the radio cone beam so it 
is not likely that the radio emission would be visible 
for this geometry. 



or cone beam. Although the force-free magnetosphere 
is a limiting case, these fits may indicate that the real 
pulsar magnetosphere solution is closer to the vacuum 
dipole in field geometry. Such light curve fitting that 
constrains the viewing geometry and radio phase lag 
will be an important tool in determining the geometry 
of the pulsar magnetosphere. 
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4. CONCLUSIONS 

We have explored the gamma-ray light curves in 
pulsar magnetosphere geometries at the two extremes, 
vacuum dipole and force-free, which both have off- 
set polar caps. In the slot gap emission model, the 
offset of the polar caps causes two major effects: a 
larger phase lag of the gamma-ray peaks with the mag- 
netic pole, and a larger ratio of peak-to-off-peak emis- 
sion level due to asymmetry in the slot gap emission. 
The sweepback of field lines of the force-free magneto- 
sphere is larger than for the vacuum dipole, producing 
a larger polar cap offset. The force-free model gamma- 
ray light curves thus show a larger radio phase lag. We 
fit simulated gamma-ray light curves using slot gap 
emission model geometry that include these effects to 
the observed Vela pulsar 30 month light curves. We 
find that the asymmetric slot gap in the vacuum dipole 
case produces a significantly better fit to the observed 
light curve since the higher accelerating electric field 
on the trailing field lines produces a lower level of 
off-peak emission. The viewing angle of ^ = 65° for 
this fit also better matches the ( from the X-ray torus 
fit. We find that the force-free model light curve fit is 
not as favorable as that of the vacuum dipole in sev- 
eral respects. While the symmetric SG in force-free 
geometry has the lowest xl^ the asymmetric SG fit 
which takes into account the offset PC is significantly 
worse. The first gamma-ray peak has a phase lag of 
~ 0.2 with the magnetic pole, which is already larger 
than the observed Vela radio phase lag of 0.14, requir- 
ing emission on only the trailing part of a radio cone. 
However, the phase of the trailing edge of a cone beam 
would be frequency dependent and the observed Vela 
phase lag is not. By contrast, the gamma-ray phase 
lag of the vacuum dipole light curves are smaller and 
are consistent with viewing the radio core beam or 
the edge of a cone beam. Furthermore, at the impact 
angle (3 = 37.5° of the force-free model fits, one would 
not expect to intercept the bright part of a radio core 
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